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In the standard model of cosmology, the present evolution of the Universe is determined 
by the presence of two components of unknown nature. One of them is referenced as "dark 
matter" to justify the fact that it behaves cosmologically like usual baryonic matter, whereas 
the other one is called "dark energy" , which is a component with a negative pressure. As 
the nature of both dark components remains unknown, it is interesting to consider other 
models. In particular, it seems that the cosmological observations can also be understood 
for a Universe which does not contain two fluids of unknown nature, but only one fluid with 
other properties. To arrive to this conclusion, we will review the observational constraints 
from supernovae of type la, cosmic microwave background, large scale structures, and the 
theoretical results of big-bang nucleosynthesis. We will try to determine constraints on this 
unifying "dark fluid" , and briefly review different possibilities to build models of dark fluid. 



In the standard model of cosmology, the total energy density of the Universe is dominated today 
by the densities of two components: the first one, called "dark matter", is generally modeled 
as a system of collisionless particles ("cold dark matter") and consequently has an attractive 
gravitational effect like usual matter. The second one, generally refered as "dark energy" or 
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"cosmological constant" can be considered as a vacuum energy with a negative pressure, which 
seems constant today. The real nature of these two components remains unknown. The reason of 
such a distinction between both dark components is mainly historical. Indeed, the cosmological 
constant was introduced by Einstein in order to justify the existence of a static Universe. Because 
this constant can appear naturally in the Einstein equations as a new fundamental constant, 
introducing it was not problematic. Moreover, no distinction was made between the dark matter 
and the usual baryonic matter. Therefore, at this epoch, the standard model of cosmology was 
not so different from the present one, with however an interesting conceptual difference: the two 
components present in the cosmological models were at that moment not considered of unknown 
nature, because, firstly no dark matter was needed, and secondly the cosmological constant could 
be considered as a fundamental constant of Nature. Thus, the Universe was finally not so strange 
at that time, and the distinction between cosmological constant and matter still holds. Later on 
important problems have appeared. First, big-bang nucleosynthesis ha s shown that the bary ons 
represent only a very small part of the matter density of the Universe Buries &: Tvler 1998l |. so 



that dark matter has to exist and to represent a large fraction of the total density of the Universe. 
Moreover, this dark matter has to be non-baryonic, and yet the standard model of particle physics 
does not provide good candidates for such a kind of matter. No particle of dark matter have 
been detected yet, and the nature of the dark matter remains uncertain. Usual candidates for 
dark matter are based on new particle physics theories, like supersymmetry, but as these new 
theories have not be en verified yet and as som e simulations based on cold dark matter seem to give 
problematic results Kazantzidis et al. 20041 ] . this answer is not definitive yet. 
A second problem arises from the cosmological constant: if it has always been constant, it means 
that its density was very small in the early Universe in comparison to the density of the other 
fluids, and appears to be dominant today only "by chance". This is the coincidence problem. 
To solve this question of coincidence, one prefers to consider that this energy is not constant, 
and the fluid which was called "cosmological constant" is now referred as "dark energy". The 
nature of this dark energy is also unknown and even stranger than usual fluids, as its pressure is 
negative today. Usual models for such a fluid are the quintessence models, which consider that the 
behavi or of the fluid can be e xplained with a real scalar field associated to a potential dominating 



today 



Peebles &: Ratra 1988l |. However, this potential is unknown, and recent observations of 



the supernova of type la seem to indicate that the ratio of pressure over density is near to -1, 
corresponding to a real cosmological constant, and can even be less than -1, what then rules out 
the involvement of usual real scalar fields and calls for even more exotic explanations. 
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To summarize, the present standard model of cosmology assumes the existence of two components 
of completely unknown nature, and the general approach to this assumption is to find constraints 
to better understand their behaviors. The usual method to find the best values of the different 
parameters of the model is to try to predict observations and to adjust the parameters to improve 
the accuracy of the predictions. In this paper, I will consider a model in which the dark matter 
and the dark energy are in fact different aspects of a same fluid, that I will call "dark fluid". 
Such an idea is simpler because one replaces two components of unknown nature by only one. 
In the following I will consider constraints from the observations of supernovae of type la, cosmic 
microwave background (CMB), large scale structures, and the theoretical predictions from big-bang 
nucleosynthesis (BBN), and I will show then that they cannot distinguish between a model based 
on two components and a model using only one. I will also use the observations to set constraints 
on the dark fluid. In a last paragraph, I will provide some ideas on what could be the nature of 
the dark fluid. 



I. BASIS OF THE MODEL AND DEFINITIONS 



If one considers a Fricdmann-Lemaitre Universe with different fluids: photons, neutrinos, baryons 
and a dark fluid, the Priedmann equations take the form: 

. a J 3 a 3 

where P and p denote the total pressure and the total density in the Universe, respectively. It will 
be assumed in the whole article that the value of a today is 1. For the dark fluid model, pressure 
and density can be expanded as: 

P = Pr + Po , 

P = Pr+Pb + PD , (2) 

where r denotes the radiation (i.e. photons + neutrinos), b the baryonic matter and D the dark 
fluid. 

In the standard model of cosmology, the equivalent expressions are: 

P — Pr ~f~ P(f ? 

P = Pr + Pb + Pdm + Pip , (3) 

where dm denotes the dark matter and the dark energy. 

For each fluid, one can write the conservation of the energy— momentum tensor in a homogeneous 
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and isotropic spacetime, which reads: 

^(Pfluid«^) = -^fluidV • 
This equation is equivalent to the second Friedmann equation. As usual, one can define a cosmo- 
logical parameter corresponding to each fluid by: 

oO _ %uid {c\ 
^%uid - ' 

where /?o is the critical density today. A cosmological parameter containing the curvature term can 
also be written: 

k 

so that one gets finally a simple equation: 

QK + ^r + ^b + ^D = l ■ (7) 

For a flat Universe, one has ^Ik = 0, and the Friedmann equations are remarkably simplified. 
Finally, one can deflne the equation of state of the fluid as: 

and one knows that for baryonic matter oji^ = 0, for radiation uJr = 1/3 and for a real cosmological 
constant uj(f) = —1. For simplicity reasons, it is considered in the following that this fluid is per- 
fect, i.e. the entropy variations and the shear stress can be ignored. For other assumptions, it is 
necessary to specify a dark fluid model. 

A basis for the dark fluid model is now defined, and a comparison with the observations of super- 
novae of type la is given in the following. 

II. CONSTRAINTS FROM THE SUPERNOVAE OF TYPE lA 

Cosmological constraints from the supernovae of type la are based on the joint observations of the 
redshift z and of the apparent luminosity / of a large number of supernova. Supernovae of type la 
are often considered as standard candles, i.e. the absolute luminosity L is approximately the same 
for every supernova (it is not completely true and recent studies corre ct the value of the absolute 



luminosity to reflect the deviations from the standard candle behavior 
it is possible to determine for each supernova the luminosity distance 



Tonrv et al. 200311 'l. so that 



L 
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This luminosity distance depends on the reddening induced by the expansion of the Universe, and 
thus can reveal the presence of the cosmological components, through the equation: 



diiz) 



1 + z 



S- 



dz' 



oo I 



(10) 



where z = a ^ — 1 denotes the redshift, S is given by 

sin X if /c > , 



S{x) 



and F is defined as 



F{z) 



X if /c = , 
sinh X if /c < , 



(11) 



(12) 



One can see that F is directly related to the first Friedmann equation through the term dz/dt = 
—a~^da/dt. If there is only one component ~ replacing the two dark components - it shall have 
the same influence on the luminosity distance - and then on the expansion of the Universe - as 
dark matter and dark energy would have. Through the Friedmann equations, it seems clear that if 



PD 

Pd 



Pdm ~)~ Pip ) 



(13) 



the dark fluid would provide the same effect on the expansion of the Universe as the two compo- 
nents. 

Observations on the supernovse of type la enable to gi ve constraints on the dark component den- 



Eiess et al. 2fln4l |. From these constraints, 



sities and on the dark energy behavior at low redshift 
it should be possible to characterize the dark fluid at low redshift. At first, the cosmological pa- 
rameter corresponding to the dark fluid can be written in function of those related to dark matter 
and to dark energy: 



D 



(14) 



At low redshift, one can consider that the equation of state for the dark energy is, at first order 
in z: 



and that the equation of state for the dark fluid can have the same form: 



(15) 



(16) 
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The observations have given constraints on the values of cj^ and lo^, and one would like to deduce 
from them constraints on co^ and coj-,. The equation of state of the dark fluid writes: 



COD 



Pd 



Po 



(17) 



Pd Pdm + P^p Pdm + Pip 

At first order in z, the dark matter density evolves like pdm = P^m'^"^ ~ Pdm(l + '^^) ■ would be 
interesting to know the behavior of p^p. Let us assume that, at first order: 



Pip = P% + pI,z ■ 

The equation of conservation of the energy-momentum tensor for each fluid satisfies: 



(18) 



di 



(P^a^) 



At first order in z, this equation becomes: 

d 



^^(p0+z(p^-3p°)) = -(a;>°) 



dt 



o„0\d{l - 3z) 



(19) 



dt 



(20) 



so that the density of dark energy reads: 



3p°(l + a.°) . 



(21) 



Then, the relation between the ratio pressure/density for the dark fluid becomes: 

, 1 p°(l + 3(l + .;°)z) 

- + "^'^L(l + 3.)+P°(l + 3(l + <).) ' 

and one can determine the value of the two first terms of the expansion: 



(22) 



+ 



(23) 



The favored values for the cosmological parameters of the usual standard model from the 



supernova of type l a 



Tegmark et al. 2004 1 are: 



Riess et al. 2004| combined with the results of other observations 



K 



h= 0.70 ±0.04 



-0.012 ±0.022 



= 0.049 ±0.012 



-1.02 ±0.19 



= 0.25 ±0.04 
0.712 ±0.044 
0.6 ±0.5 . 



(24) 
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From these values, one can calculate the parameters of the dark fluid: 

n% = 0.962 ±0.084 

uj% = -0.76 ±0.25 (25) 
= 1.0 ±0.6 . 

These values are of course not completely representative of the dark fluid model, because they 
come from data analyses based on the usual standard model. Nevertheless, one can use them as 
test-parameters at low redshift. 

Recent supernova observations tend to show that the dark energy has a negative pressure. More- 
over, ojip < —1 is not at all excluded, and in that case the dark energy cannot be explained 



Caldwell et al. 2003l | for a possible answer 



anymore thanks to the usual models (see for example 
to this problem). From the precedent constraints, one can see that this difficulty vanishes with a 
dark fluid. Hence, the pressure of the fluid has to be negative today at cosmological scales, and 
seems to increase strongly with the redshift. As Wjj > —1, it seems possible to model the dark 
fluid with a scalar field. 

The study of supernovae provided us properties of the equation of state of our dark fluid at low 
redshift independently from the specification of a dark fluid model. We will now try to extract 
constraints from the information concerning large scale structures. 

III. LARGE SCALE STRUCTURES 

We will not consider here a complete scenario of structure formation, which would require the 
specification of a precise model of dark fluid. Nevertheless, one can study the necessary conditions 
for the fluid parameters to enable the perturbations to grow and to give birth to large scale 
structures. 

Let us consider the case where the equation of state of our fluid does not change during the growth 
of perturbations, and, to simplify, that the entropy perturbations can be ignored and that the 
Jeans length is smaller than any other considered scale. One can define the local density contrast 
of the dark fluid as: 

o{x,t) = — 1 , (26) 

PD(t) 

where pD {x, t) is the local value of the density, and 'PD{t) is the mean background density, i.e. the 
apparent cosmological density. In the fluid approximation, one can write the evolution equation of 
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the local density contrast Gaztanaga &; Lobo 200ll |: 



I dF5 ( H \ 1 dd 2 
7 

4: + 3uD 1 / 1 \2 1 



TTo bl + o (1 + ^d){1 + 3^z5)f7z5<5^ . (27) 



3{l + lod) 1 + S \ l + 5 J H'^ydtJ 2 
To solve this equation, one can define a new variable reflecting the expansion: 

?? = lna , (28) 

so that, if one assumes that the dark fluid is completely dominant at the time of growth of pertur- 
bations (in that case, the Friedmann equations reveal that H/H^ = —3(1 + a;/))/2), equation (|27|1 
becomes: 

(29) 

Because the coefficients of the above equation are time-dependant only, one can separate the spatial 
and temporal parts so that 

Si{x,t)=5o{S)D{t) , (30) 

where D is called the "linear growth factor" . In the linear approximation, where 6 is small, equation 
becomes: 

d'^D l-3ujDdD 2, 

■W + - 3(1 + + ^o)D = . (31) 

Its solutions take the form D = Dia"^ + D2a"^, with Di and D2 being two integration constants, 
and 

ai = 1 + 3u;d , (32) 
a2 = -^(I+cod) . (33) 

Thus, in the case of a dominant dark fluid, we have only a growing mode if uju > —1/3 or if 
uj£) < —1. One can however note that the last inequality seems very difficult to achieve with 
standard model for dark matter or dark energy models. 

Let us now consider the observations of the cosmic microwave background to get constraints at 
earlier times. 
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IV. COSMIC MICROWAVE BACKGROUND 

A power spectrum of temperatu re fluctuations can b e deduced from the observations of the cosmic 
microwave background (CMB) 



Spergel et a 



hard work, and a program like CMBFAST 



Predicting this p ower spectrum requires a 



Seliak k. Zaldarriaea 1996l | is able to produce it for 



the cosmological standard model. In our case, we will consider only the position of the peaks to 
constrain the parameters of the dark fluid, and we will make some assumptions on the dark fluid 
properties. 

First, one should note that at high redshift, in the standard model the density of dark energy is 
nearly negligible in comparison to that of the dark matter. As the usual model seems to be able 
to correctly reproduce the fluctuations of the CMB, one can assume that our dark fluid should not 
behave very differently from the superposition dark matter/dark energy, and so should behave at 
the moment of recombination nearly like matter. Therefore, one can write the density of our fluid 
as a sum of a matter-like term (m) and of another term of unknown behavior (o): 

, n ^ -3 



PD = pDm\ — ] +PDo ■ (34) 

One can note that this equation gives no constraint on the behavior of the dark fluid, as the second 
term is not restricted to any behavior yet. 

We do not want to specify a model of dark fluid and we would like to be as general as possible. 
Nevertheless, we will consider for simplicity only the background properties of the dark fluid, and 
we will not try to reproduce the whole power spectrum, but only consider the position of its peaks 
without trying to find their amplitude. The conformal time is defined by: 

dt a-^{t) . (35) 



The spacing between the peaks is then given, to a good approximation, by 



Hu fc Sutrivama 199,1| : 



A/«7r^^-^ , (36) 

CsTls 

where Cs is the average sound speed before last scattering, and tq and tis the conformal time today 
and at last scattering. This average sound speed reads: 

"1/2 



where p}, is the density of baryonic matter and pr is the density of relativistic fiuids (radiation and 
neutrinos). 
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Let us consider that the Universe is flat so that the Priedmann equations are simpHfied. 
In this case, the first Priedmann equation can be written: 

2 SttG . , , , 

H = {Pb + Pr + PDm + PDo) ■ (38) 

Using the evolution equation of the different densities, this becomes: 

= Hi (n^a-^ + n'^a-' + n^^,^ (^) + ^-^pno • (39) 

The precedent equation cannot be solved if the form of pDo is not given. In our case, it is possible 
to assume that the fraction 

does not vary too rapidly before the moment of last scattering (denoted Is), so that an effective 
average can be defined: 



^Do 



^ r^s' r ^Do{r)dT . (41) 

^0 



In the following, we will only consider an approximate and effective density: 



SttG 

Let us then replace this density in the Friedmann equation: 



PDo^H^^d^o- (42) 



H\l-n'^,) = Hl[{nl + n'i,^al)a-' + n%-^) . (43) 

This time, provided one fixes the values of the different cosmological parameters and knowing the 
initial conditions, this equation can be solved. While replacing usual time by conformal time, the 
Priedmann equation becomes: 

{t) ' = - (^^° + ^Dmai)a{r) + ^f) . (44) 

The value of the conformal time at the moment of last scattering is given by: 



One can use the same method to evaluate the conformal time today. The Priedmann equation 
reads, after last scattering: 



{^y = Hi[nUr) + n'r+airr ^) • (46) 
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Let us make now the further assumption that the dark fluid has an approximate equation of state: 

where is an effective value of the dark fluid density, such that 

and ud is the average value oi ud over the conformal time, weighted by: 

to reflect the fact that the equation of state of our fluid should be more significant when its density 
contributes more heavily to the total density of the Universe, ojd is then given by: 

/■To 

UD = . (50) 

/ ^D{T)dT 

Jo 

Defining the effective cosmological parameter O^j = ^P%/{^t^GHq), the Priedmann equation be- 
comes: 

(^) ' = Hi {n1a{r) + n'^r + n%air)('-^-^) . (51) 
One can then integrate the equation, and show that: 

to = 2Ho^F{Ud) , (52) 

with 

F{uJd) = ^£ da(nla + n°, + Cl%a^^-'^'^^y^^^ . (53) 

There is no analytical integration for this function, but in a few cases. In particular, one has for a 
flat Universe: 

^(0) = n^^lci% {^^' + ^'r + ^%- V^) , (54) 

and 

F{-m = lin%)-'/'lnl ^ ' , y . (55) 

2 \9P, + 2,J^^] 
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The case lod = —1/3 looks much more probable than ujd = when one considers the value of lo^ 
which was obtained from the supernova data. One finally gets the spacing between peaks: 



F(uJd) 



'■Do 



The sound velocity c<j is then given by: 



17° 



Is . 



(56) 



-1/2 



(57) 



This equation can be integrated analytically, and one finally gets: 

1/2 f I — I — N -1 



no 



-c -if^ 



One can notice that Cs does not depend on r^/jo- The approximate value of ais can be taken 



from 



Hu et al. 2nnil| : 



aj-,^ ^ 1008(1 + 0.00124(0>2)-o.74)(;l ^ ci(Sl^„a? 



3 \C2^ 



(59) 



where 



ci = 0.0783(Q^/i^)-"-^^(l + 39.5(0^/1^) 
C2 = 0.56(1 + 21. l(l^|]/i2)i-28)-i . 



0i,2n0.76\-1 



(60) 



Unfortunately, this is an implicit equation, and we cannot find ais analytically. 
Finally, one has a direct dependence between Al and the parameters ^^^m, ^^o' ^^'^ ^D- Even 
if this is only an approximate formula, it gives the possibility to visualize directly the effect of the 
dark fluid on the position of the peaks, provided that its behavior does not differ much from the 

— Is 

requirements of the approximations on $7^)^ and cj/j. 

However, the calculated Al cannot be directly related to the observed spacing between peaks, as 
shifts of peaks can be induced by other effects. In particular, the location of the i-th peaks can be 
approximated by: 



= Al(m 



Al(m 



(61) 



where <j) is the shift of the first peak, corresponding to an overall shift, and the 5(i)i is the specifi c 



shift of the i-th peak. We have fortunately access to the fitting formulae of 



Doran fc Lillev 2002l| : 
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Overall phase shift 



1.466 - 0.466n) airf^ + 0.291J7£,„ 



where ai and 02 are given by: 



ai = 0.286 + 0.626 
a2 = 0.1786 -6.308Slfe/i2 + 174.9 (l^fo/i^)^- 1168 



n is the spectral index and ri^ is defined by: 

Priais 



ris 



PDm{ais) ^n, 



af 



Dm"'ls 



Relative shift of second peak 6(j)2 : 



Co - ciHs - C2r,„''^ + 0.05 (n - 1) , 



with 



Co = -0.1 + (0.213 - 0.12317^ exp {- (52 - 63.617^o) , 
ci = 0.063 exp I -3500 (Jlfo/i^)^! +0.015 , 
C2 = 6 X 10'^ + 0.137 (Ofc/i^ -0.07)^ , 



C3 



0.8 + 2.3S7^o + (70 - 12617^ J^fe/i^ . 



Relative shift of third peak 



63 = 10 - di rfj^ +0.08 (n - 1) 



with 



di = 9.97+ (3.3 -317^„) l^b/i^ , 



(62) 



(63) 



(64) 



(65) 



(66) 



(67) 



(68) 



d2 = 0.0016 - 0.0067n^, + (0.196 - 0.22^^,) Q.h^ + ^^'^^ + ^'IJ^?"^ 

One can now compare our results to the da ta. The WMAP experiment provides the precise location 
of the two first peaks 



Suergel et al. 2003l| : 



220.1 + 0.8 , 



Ip^ = 546 + 10 , 



(69) 
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' 'Dm 


"Do 




Ud 


h 


h 


^3 


0.600 


0.125 


0.95 


-0.10 


220.3 


544.8 


817.2 


0.400 


0.120 


0.95 


-0.15 


220.0 


546.1 


830.6 


0.270 


0.080 


0.95 


-0.20 


219.4 


544.8 


835.4 


0.200 


0.010 


0.50 


-0.13 


220.2 


544.5 


835.9 


0.250 


0.035 


0.50 


-0.10 


220.7 


544.1 


833.5 


0.018 


0.013 


0.25 


-0.07 


219.6 


543.4 


837.4 


0.017 


0.001 


0.10 


0.17 


219.5 


543.1 


837.6 



TABLE I: Position of the peaks for h 0.7, n ^ 1, n° ^ 0.049 and Q'^. = 9.89 x 10"^ in function of the 
parameters of the dark fluid. 



and BOOMERanG gives the position of the third peak Bernardis et al. 20n2l |: 



Ip^ = 825 ± 13 . (70) 

To evaluate roughly the value of the parameters of the fluid, one can fix the other parameters as 
follows: 

n = 1 , 

h = 0.70 , (71) 
= 0.049 , 
Q° = 9.89 X 10"^ . 

For these values, one finds in Table 1 the resulting positions of the peaks in function of parameters 
of the dark fluid. 

One can note that a large range of values is possible. It is not so strange because we have many 
parameters for our fluid. A more complete analysis is not needed here, because the strongest 
constraints would come from the specification of a model. Without specifying a model, Table 1 
shows that the values of the parameters of the dark fluid are not stringently constrained. One can 
nevertheless see that for large values of Cl^, the permitted values of tZJ^ are negative, and hence 
in that case one can assume that our fluid behaves today like a cosmological constant whereas it 
could have behaved mainly like matter at last scattering. For small values of Cl^, Ud is positive, 
so that the density of dark fluid should decrease more rapidly than a matter density after last 
scattering. In this case, fi^o is very small, and the fluid should have behaved like matter before 
and around last scattering. Small values of Cl^ look therefore unrealistic, because as uJd is then 



15 



positive, unless our dark fluid has an oscillating density, i^Do should be much larger and certainly 
dominant. 

— Is 

The value of i^/jo, which can be as much as 0.1, also shows that before recombination, the fluid 
may have behaved differently from matter, and perhaps like radiation. One can also note that low 
values of corresponds to low values of ^2^^, and then in that case the fluid mostly behaves like 
matter. 



If one combines these results with the results from the CMB, it seems that the constraints 
become — l/3<u;z)<0to enable the perturbations to grow. This shows in fact that during the 
growth of the perturbations the behavior of the dark fluid should not be too different from that 
of matter. This result also confirms that the value of the effective which appears in equation 
(|51jl cannot be too small (it has to be at least larger than 0.2). 

If one wants to perform a much more precise study of a specified model, it would be interesting 
to sim ulate the whole process of structu re formation, and to compare the results with surveys like 



SDSS 



Tegmark et ab 2004 Il| . or 2dF 



Percival et al. 200 ll |. We will not study here further the 



CMB power spectrum, as other features seem more model-dependant, and we want to consider 
here the general case. We will now consider the results of the big-bang nucleosynthesis and their 
influence on the establishment of a dark fluid model. 



V. BIG-BANG NUCLEOSYNTHESIS 



Recent analyses of the big-bang nucleosynthesis (BEN) |Coc et al. 2004l | indicate a discrepancy 
between the value of the baryonic density calculated from the observed Li and ^He abundances, and 
the one calculated with the observations of deuterium. Some explanations can be found. Problems 
could have appeared in the measurement of Li and "^He abundances, or the Li on the stellar surface 
could be altered during stellar evolution, or we have no accurate knowledge of the reaction rates 
related to ''Be destruction, or the expansion rate during BBN could have been modified through 
an accelerated cosmological expansion. Thus, two possibilities can be considered for the equation 
of state of the dark fluid. First, if it is correct to consider a Universe dominated by radiation at 
BBN time, the main constraint is that the dark fluid density should be small in comparison to the 
radiation density; otherwise Friedmann equations indicates that the expansion rate of the Universe 
would be different from the one in the standard BBN, changing then the temperature evolution 
rate and so the abundance of the elements. It means that, if one assumes that the dark fluid 
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behavior does not change violently during BBN, the equation of state of this fluid around the time 
of BBN has to be C(;£)(BBN) < 1/3, or that its density was completely negligible before BBN. In 
the case of a real radiative behavior W£)(BBN) = 1/3, the dark fluid behaves like extra-families of 
neutrinos, and its density can be constrained. The effective extra-neutrinos number at the BBN 
time is defined by: 

^^e.(BBN) . « . (72, 

where py is the standard density of a single relativistic ne utrino species. The usual bound on 



Buries &: Tvler 1998l |. which corresponds in our case for a 



the number of neutrinos is IS.N^^ < 1 
temperature around 1 MeV (a^^ ~ 3 — 4 x 10^) to: 

7 / 4 X 4/3 ^2 

pz?(BBN) < - (^-j Y5^BBN ~ ^ x 10-2(MeV)^ . (73) 
This limit is only valid in the case a;£)(BBN) = 1/3. 

If the abundance of the elements is as observed, a modification of the expansion rate could 
provide, as presented in [Salati 20031 ]. a correction to the predicted values. If our fluid is the 
dominant component at the time of the BBN, it can have a big influence on the expansion rate 
of the Universe. Evaluating the density of the dark fiuid and its evolution during the time of the 
BBN so that the observations are retrieved would require further studies that I will not develop here. 



We have seen that a dark fiuid may be compatible with the cosmological observations, and 
could be an interesting approach to the ambivalence of dark energy and dark matter. Let us now 
consider possible paths to model a dark fiuid. 



VI. MODELS OF DARK FLUID 



In the literature, only few fiuids behaving like a dark fiuid are considered. Different ways to model 
the dark fluid are possible, and I will consider here in particular two of them: the generalized 
Chaplygin gas, based on D-brane theories, and another one using scalar fields. 



A. Generalized Chaplygin Gas 



The generalized 



Chaplygin gas (GCG) is an exotic fiuid derived from D-brane theories 



Bento et al. 2002l |. It can be described by an equation of state: 



ch 



A_ 



(74) 
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where a is a constant, < a < 1, and A is another positive constant. This equation of state 
corresponds to a density evolving hke: 



Pch — Pch ( H gTY 



(75) 



a3(l+a) 

where Ag = ^/(/j^!^)^-'^^") and pj?^ is the Chaplygin gas density today. Such a behavior could be 
interesting in order to model the dark fluid, because for high values of a this density is mainly 
constant, and for low values of a it evolves like matter. This behavior has to be compared with 
the observations. 

For the comparison with the data of supernovae of type la, let us derive the equation of state of 
the GCG at low redshift: 

Pch 



"""^ Pch ^. + (l-^)(l + z)3(l+° 

SO that one can deduce: 



-^[l-3(l-^)(l + a)z] , 



(76) 



^'ch 



-A.. 



ujlh = 3^(l-^,)(l + a) 



(77) 



to be compared with the constraints (|25j) . and one finally gets: 

As = 0.76 ±0.25 , 
a = 0.8 ±2.4 . 

We have of course no constraint on a. The analysis of CMB has been done in 



(78) 



Bento et al. 200311 . 



and the results for h = 0.7 and n = 1, combined with the ones from the supernova of type la are: 



As = 0.75 ±0.13 
a = 0.4 ±0.2 . 



(79) 



If one considers now the BEN, at this time the GCG behaves like matter, so that it is compatible 
with the standard BEN scenario. The large scale structure formation has been studied in a case 
where the Chaplygin ga s adds only a background density to a Universe containing cold dark matter 



Multamaki et al. 2004l |. and such a scenario seems then possible. 



So, the GCG seems to be in agreement with the cosmological observations. A further analysis 
is still needed, in particular concerning the growth of structures with a dominant Chaplygin gas 
density, or the local behavior of such a fluid. 
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B. Scalar Fields 

One can also consider the idea that the dark fluid could be explained thanks to a scalar 
ield. Indeed, scalar f i elds are very useful in expla ining the behavior of the dark energy today 



Eatra fc Peebles 2nnr 



can he 
scales 



lave like matter on local scales 



Arbev et al. 20021 



Hebecker fc Wetter ich 2001 1, and recent analy ses have shown that they 



Arbev et al. 200.'^ 
Guzman fc Urena-Lopez 200,j| 



Kiselev 2005l | as well as on cosmological 



Let us therefore consider a real scalar field associated with a Lagrangian density 



d^ifduLp - V {lp) . 



(80) 



Its pressure and its density on cosmological scales are given by 



1 .2 



1 



v3^ + y(^) . 



(81) 



So, the pressure is negative if the potential dominates, and negligible if the potential equilibrates 
the kinetic term. Thus, a scalar field can be a good candidate for the dark fluid if it respects in 
particular the following constraints: 

- its density at the time of the BBN decreases at least as fast as the density of radiation, 

- its density from the time of last scattering to the time of structure formation evolves nearly 
like matter, and so ^v?^ « V {(p), 

- after the growth of perturbations, because the scalar field is dominating the Universe, its 
potential does not equilibrate the kinetic term anymore and begins to dominate; thus it will 
behave like a cosmological constant in the future. 



The main parameter of such a model is the same as that of quintessence models: the po- 
tential. When one considers quintessence models with real scalar field, one looks for potentials 
which provide a cosmological constant-like behavior today, and dec reasing potentials s eem to be 



favored. If one considers now complex scalar fields, it was shown in [Arbev et al. 2002l | that such 
a field can behave like cosmological matter when its potential has a dominant m^|(/)p term, which 
corresponds to an increasing term in the potential. Thus, a way to find a "good" potential would 
be to consider a superposition of a decreasing potential which would begin to dominate today, 
and of the increasing quadratic term which has to dominate at least until structure formation and 
can nevertheless lead to an attractive effect on local scales today. A more detailed study of the 
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possible potentials will be presented in a future publication. 



One could also hope that a scalar field might explain the excess of gravity on local scales. 
To do that, one can imagine a Universe filled with a scalar field. In the part - and time - of the 
Universe where the density of baryonic matter is high, the scalar field would, through gravitational 
interaction, have a large kinetic term which could even equilibrate the potential, so that one has 
an attractive net force o n local scales (easie r to achieve with a complex scalar field associated to 



an internal rotation, see 



Arbev et al. 2003l |). whereas in the parts where no baryons are present. 



the field would not vary much, and the potential dominates, providing repulsion. Thus, on local 
scales, where the baryon density is high, the field behaves like matter. Where the baryonic density 
is small, i.e. away from galaxies and clusters, the gravitational interaction is not strong enough 
to increase the kinetic term of the scalar field, so that the potential dominates, and one can 
then observe the effects of a negative pressure. In that case, the scalar field will have a negative 
pressure on cosmological scales, providing a locally negative pressure in average. In the past, 
baryons were uniformly dense, so that the kinetic term was large everywhere, and one could have 
then a uniform matter behavior under these conditions. In that way, the local behavior can be in 
agreement with the cosmological one, and a complete cosmological scenario can be built. Such a 
scenario has of course to be studied further. 



VII. CONCLUSION AND PERSPECTIVES 



Astrophysical and cosmological observations are usually interpreted in terms of dark matter and 
dark energy. We have seen here that they can also be analyzed differently. Thus, it is possible 
to develop a model of dark fluid, which could advantageously replace a model containing in fact 
two dark components. Of course, hard work and studies are required to test completely the dark 
fluid hypothesis. Nevertheless today, as it seems difficult to find a model for dark energy and 
as problems concerning cold dark matter remain, it is worthwhile to investigate different ideas 
such as an unification of dark energy and dark matter - that finally does not seem stranger than 
trying to determine the nature of two components at the same time - which could be achieved in 
particular thanks to D-brane theories (in particular through the Chaplygin gas), or thanks to the 
so-useful scalar fields. Of course, other models may also account for dark fluid. 



An important question remains, how to interpret the dark matter problem on local scales 
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and could the dark fluid account for the excess of gravity inside local structures? I will only pro- 
vide here a qualitative analysis of whether a fluid with a negative pressure on cosmological scale can 
have an attractive effect on loc al scale, such as it i s observed in galaxies ffor example, with the r o- 



tation curves of spiral galaxies 



Persic et al. 1986 , 



Gentile et al. 2004 . 



Carignan Purton 1998| |'). 



Let us consider only the quasi-Newtonian limit of general relativity. In that case, devia- 
tions from the Minkowski metric rj^u = diag(l, —1, —1, —1) are accounted for by the perturbation 
/i^j/. In the harmonic coordinate gauge, it satisfles the condition: 

dah'^, - \d^hl = . (82) 

One can show that the perturbation is related to the source tensor: 

S/j,!/ = T^u ~ 2 dfiiy T a (83) 



V(f) = -4G / ^If^dV . (84) 



through the integral 



If the energy-momentum tensor is written as: 

rpi^u = (p + pp/'U" - Pg'"' , (85) 

one can show that, at leading order in h^^^, 

S^''' = (P + p)U^'U'' -^q^^ip- P) . (86) 

The gravitational potential is in fact <I> = /ioo/2, so that: 

$(f) = -2Gjp^^d^r' . (87) 

For a fluid at rest, = (1, 0, 0, 0), with an equation of state P = cop, one has: 

Soo = + 3a;) , (88) 

and consequently this fluid has an attractive effect only if u; > —1/3. From the study of super- 
novae it seems that our dark fluid is not in that state today, so that its effects are mainly repulsive. 
Nevertheless, it is possible that the dark fluid has a different behavior on cosmological and on local 
scales. Indeed, the density and pressure of the fluid on cosmological scale are spacial averages of 
the local density and pressure, and one can assume that: 

p{t,r) = p^^^^^{t)+6p{t,r) , 

P(t,f) = P^°^™°(i) +<5P(t,f) , (89) 
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where p^osmo g^j^^j pcosmo ^j.^ ^^le cosmological density and pressure, with the spatial averages: 

< 6p {t, r) >=< 6P {t, r) >= . (90) 

The density of dark fluid on cosmological scales today is of the order of the critical density, i.e. 
p^! ~ 9 X 10~^^g.cm^^. One can compare it to the estimated matter den sity in the Milky Way at 
the radius of the Sun pSun ^ 5 X IQ-^^g.cm-^ |oiling fc Merrifield 200^ . Hence, even if the dark 



fluid's local density would represent 1% of this total matter local density only, its value would be 
much higher than the cosmological densities today. Therefore, on local scales, one can assume that 
5p {t, r) » pCosmo (-^^ ^^m^ write: 

So^^\{5p + ^5P) . (91) 

To have a net attraction, we get finally the same kind of constraint as before, 6p > —36P, but this 
time we do not have to use the cosmological constraints, because the local behavior of the dark 
fluid can be very different from the cosmological one. Moreover, if no model is specified, one can 
still hope that 5P could be negligible on local scales, so that the local behavior of the dark fluid is 
matter-like, and that the usual Newtonian equation can be retrieved: 

$(^-^ = -G [ -M^d^r' . (92) 
J \r' — r \ 

This local behavior will of course have to be verified quantitatively for each dark fluid model, but 
nevertheless gives hope for a unified explanation on any scale. In particular, considering scalar 
fields, this qualitative analysis tends to show that it would be interesting to try to find a potential 
which gives today a negative pressure on cosmological scale, but which also gives a matter behavior 
in local structures, i.e. where the density of baryons is high. 

To conclude, the dark fluid appears as an interesting possibility to explain the observations. As 
the properties of the dark fluid are different from dark matter and dark energy, models of dark 
fluid are worth to be studied, and we can now use many precise observations as strong constraints 
on such models. 



22 



Acknowledgements 

I would like to thank Farvah Mahmoudi, Helene Courtois, Julien Devriendt, Thierry Sousbie and 
Wolfgang Hillebrandt for their comments and for useful discussions. 



[Arbey et al. 2002] A. Arbey, J. Lesgourgucs & P. Salati, Phys. Rev. D 65 (2002), 083514. 
[Arbcy ct al. 2003] A. Arbey, J. Lesgourgues & P. Salati, Piiys. Rev. D 64 (2001), 123528; Piiys. Rev. D 68 
(2003), 023511. 

[Bcnto et al. 2002] M.C. Bento, O. Bertolami & A.A. Sen, Phys. Rev. D 66 (2002), 043507. 
[Bento et al. 2003] M.C. Bento, O. Bertolami & A.A. Sen, Phys. Lett. B 575 (2003), 172. 
[Bernardis et al. 2002] P. Bernardis et al., Astroph. J. 564 (2002), 559. 

[Buries & Tyler 1998] S. Buries and D. Tyler, Astrophys. J. 499 (1998), 699; Astrophys. J. 507 (1998), 732. 
[Buries & Tyler 1998] S. Buries & D. Tyler, Astrophys. J. 499 (1998), 699; Astrophys. J. 507 (1998), 732. 
[Caldweh et al. 2003] R.R. CaldweU, M. Kamionkowski & N.N. Weinberg, Phys. Rev. Lett. 91 (2003), 
071301. 

[Carignan & Purton 1998] C. Carignan & C. Purton, Astrophys. J. 506 (1998), 125. 
[Coc ct al. 2004] A. Coc et al., Astrophys. J. 600 (2004), 544-552. 

[Doran & Lilley 2002] M. Doran & M. Lilley, Mon. Not. Roy. Astron. Soc. 330 (2002), 965. 

[Gaztanaga & Lobo 2001] E. Gaztafiaga & J. A. Lobo, Astrophys. J. 548 (2001), 47. 

[Gentile et al. 2004] G. Gentile et al., Mon. Not. Roy. Astron. Soc. 351 (2004) 903. 

[Guzman & Urcna-Lopez 2003] F.S. Guzman, L.A. Urena-Lopez, Phys. Rev. D68 (2003), 024023. 

[Hebecker & Wctterich 2001] A. Hebcckcr & Wettcrich, Phys. Lett. B 497 (2001), 281. 

[Hu et al. 2001] W. Hu, M. Fukugita, M. Zaldarriaga & M. Tegmark, Astrophys. J. 549 (2001), 669. 

[Hu & Sugiyama 1995] W. Hu & N. Sugiyama, Astrophys. J. 444 (1995), 489. 

[Kazantzidis ct al. 2004] S. Kazantzidis ct al., Astrophys.J.608 (2004), 663. 

[Kiselev 2005] V.V. Kiselcv, Class. Quant. Grav. 22 (2005), 541. 

[Muhamaki et al. 2004] T. Muhamaki, M. Mancra & E.Gaztafiaga, Phys. Rev. D 69 (2004), 023004. 
[Olhng & Merrihcld 2001] R.R Oiling & M.R. Merrifield (Nottingham U.), Mon. Not. Roy. Astron. Soc. 
326 (2001), 164. 

[Peebles & Ratra 1988] P.J.E. Peebles & B. Ratra, Astrophys. J. 325 (1988), L17. 
[Percival et al. 2001] W.J. Percival et al., Mon. Not. Roy. Astron. Soc. 327 (2001), 1297. 
[Persic et al. 1986] M. Persic, P. Salucci & F. Stel, Mon. Not. Roy. Astron. Soc. 281 (1986), 27. 
[Ratra & Peebles 2000] B. Ratra & P.J.E. Peebles, Phys. Rev. D 37 (2000), 3406. 
[Riess et al. 2004] A. G. Riess et al, Astrophys. J. 607 (2004), 665. 
[Salati 2003] P. Salati, Phys. Lett. B 571 (2003), 121. 



23 



[Spergel et al. 2003] D.N. Spergel et al., Astrophys. J. Suppl. 148 (2003), 175. 
[Seljak & Zaldarriaga 1996] U. Seljak & M. Zaldarriaga, Astrophys. J. 469 (1996), 437. 
[Tegmark et al. 2004 I] M. Tegmark et al., Phys. Rev. D 69 (2004), 103501. 
[Tegmark et al. 2004 II] M. Tegmark et al., Astrophys.J.606 (2004), 702. 
[Tonry et al. 2003] J.L. Tonry et al., Astrophys. J. 594 (2003), 1. 



